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The angular distributions are studied in events with at least two jets; the highest dijet mass observed 
is 5.5 TeV. All angular distributions are consistent with the predictions of the Standard Model. In a 
benchmark model of quark contact interactions, a compositeness scale below 8.1 TeV in a destruc¬ 
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CL; median expected limits are 8.9 TeV for the destructive interference scenario and 14.1 TeV for the 
constructive interference scenario. 
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Search for New Phenomena in Dijet Angular Distributions in Proton-Proton 
Collisions at y / s = 8 TeV Measured with the ATLAS Detector 

ATLAS Collaboration 
(Dated: 31 March 2015) 

A search for new phenomena in LHC proton-proton collisions at a center-of-mass energy of y/s = 

8 TeV was performed with the ATLAS detector using an integrated luminosity of 17.3 fb -1 . The 
angular distributions are studied in events with at least two jets; the highest dijet mass observed is 
5.5 TeV. All angular distributions are consistent with the predictions of the Standard Model. In a 
benchmark model of quark contact interactions, a compositeness scale below 8.1 TeV in a destructive 
interference scenario and 12.0 TeV in a constructive interference scenario is excluded at 95% CL; 
median expected limits are 8.9 TeV for the destructive interference scenario and 14.1 TeV for the 
constructive interference scenario. 


The search for an internal structure of fermions and 
new forces that might govern that structure is a ma¬ 
jor goal of modern particle physics. The most powerful 
probes are scattering experiments with large momentum 
transfer. Collisions of protons at the Large Hadron Col¬ 
lider (LHC) resulting in two energetic jets of particles 
(dijets) provide the largest momentum transfer currently 
available, and therefore the deepest probe. 

The angular distribution of jets relative to the beam 
axis in events with high dijet invariant mass ( rrijj ) pro¬ 
vides stringent tests of perturbative quantum chromo¬ 
dynamics (QCD) as well as theories of new phenomena. 
QCD calculations predict that dijet production, domi¬ 
nated by t-clrannel gluon exchange in the kinematic re¬ 
gion of proton-proton (pp) collisions at the LHC, fea¬ 
tures steeply falling rrijj distributions and angular distri¬ 
butions peaked at |cos(#*)| = 1, where 0* is the polar 
scattering angle in the two-parton center-of-mass frame. 
New phenomena, such as strong gravity mm or new 
interactions typically predict angular distributions 
which are more isotropic. 

Previous studies of dijet angular distributions, at the 
CERN SPS [7UH], the FNAL Tevatron [1 [TO] . the CERN 
LHC at y/s = 7 TeV (TTl - fTfi] . and by the CMS Collabo¬ 
ration at y/s = 8 TeV [17], have reported results consis¬ 
tent with the Standard Model (SM). This Letter reports 
on studies of dijet angular distributions in pp collisions 
at y/s = 8 TeV in data with an integrated luminosity of 
17.3 fb -1 collected with the ATLAS detector in 2012. 

A detailed description of the ATLAS detector is pub¬ 
lished elsewhere m • The detector is instrumented over 
almost the entire solid angle around the pp collision point, 
with layers of tracking detectors, calorimeters, and muon 
detectors. 

The jets are measured using a calorimeter system com¬ 
posed of different detector types covering different regions 
in rj BS] and depth. The electromagnetic calorimeter 
is composed of liquid-argon sampling calorimeters using 
lead as absorber, and is split into a barrel (|?y| < 1.475) 
and two end caps (1.375 < |?y| < 3.2). The hadronic 
calorimeter is divided into a barrel and two extended 
barrels (|p| < 1.75), and two endcaps (1.5 < \r/\ < 3.2). 
The barrel and extended barrels are sampling calorime¬ 


ters with steel as absorber and scintillator tiles as the 
active medium, while the hadronic endcaps are liquid- 
argon calorimeters with copper as absorber. In the very 
forward regions (3.1 < |ry| < 4.9) there are liquid argon 
calorimeters with copper and tungsten absorbers. 

The data are selected using a trigger that requires a 
single high-px m jet above one of eight thresholds, rang¬ 
ing from 25 GeV to 220 GeV. Due to the high rate of jets 
at lower px, only a fraction of the events from the lower 
seven thresholds are stored. 

Individual jets are reconstructed using the anti-fct 
jet clustering algorithm [5l3 ED with radius parameter 
R = 0.6. The inputs to this algorithm are clusters [22] 
of calorimeter cells with energy depositions significantly 
above the noise. Jet four-momenta are constructed by 
the vectorial addition of clusters of cells, treating each 
cluster as a four-momentum with zero mass. The jet four- 
momenta are then corrected to the jet energy scale [23] 
as a function of ?y and px for various effects, the largest 
of which are the hadronic shower response, detector ma¬ 
terial distribution, and pileup events [24j . This is done 
using a calibration scheme based on samples of simulated 
events and validated with test-beam [25| and collision 
data [22] studies. 

The rapidity of a jet is defined as 
y = \\w[{E + p.)/(E — p z )\, where E is the jet en¬ 
ergy and p z is the momentum component along the 
beam axis [T5]. The scattering angle between two jets 
can be expressed using the variable y = e^ Vl ~ V2 ^ = e 2 \ y I, 
where j/i and y 2 are the rapidities of the two jets, and 
y* = \{y i — y 2 ). The rapidity boost of the dijet system 
with respect to the center of mass of the colliding 
protons is calculated as ps = \{yi + 2 / 2 )- 

For each trigger, the event is required to have a jet 
with px sufficient to achieve a trigger efficiency greater 
than 99.5%. For the lowest (highest) threshold trigger, 
this corresponds to px > p“ ln = 47 (333) GeV. Events 
are required to have at least two jets, each with px > 
50 GeV; the dijet system, defined as the two jets with 
largest px, is required to have |p*| < 1.7, |pe| < 1-1, and 
rrijj > 600 GeV (where the rrijj requirement avoids the 
kinematic bias in the angular distributions introduced by 
the minimum px requirement). 


The detector covers the angular range \y*\ < 1.7, corre¬ 
sponding to x < 30. This interval is divided into 11 bins, 
with boundaries at Xn = e^°- 3xrl ^ for n = 0 to 10, ap¬ 
proximating the segmentation of the calorimeter in Arj. 
The data are further binned coarsely in rrijj with the 
expectation that low-m 7J - bins are dominated by QCD 
processes and that signals associated with new physics 
would be found in higher dijet invariant mass bins. The 
bin edges are chosen to optimize the expected sensitivity 
to the model of contact interactions. The highest dijet 
mass observed is 5.5 TeV. 

The SM predictions are estimated using the 
Pythia8 |26j v8.160 event generator with the AU2 [27] 
underlying-event tune and the CT10 [23] parton dis¬ 
tribution functions (PDF). The simulated events are 
propagated through a detector simulation [25j that uses 
the Geant4 [30] simulation package. Pileup conditions 
vary as a function of the instantaneous luminosity 
and are taken into account by overlaying simulated 
minimum-bias events generated with Pythia8 onto 
the hard-scattering process such that the observed 
distribution of the average number of interactions per 
bunch crossing is reproduced. The same reconstruction 
and event selection are applied to the simulated events 
and the data. 

The Pythia8 calculations are primarily to leading or¬ 
der (LO) in QCD with simulation of higher-order contri¬ 
butions included in the shower modeling. Events gener¬ 
ated by Pythia8 are reweighted using a correction factor 
calculated based on the ratio of the next-to-leading-order 
(NLO) cross-section calculation from Nlojet-| —b [3TH33] 
v4.1.2 to the LO+shower calculation from Pythia8: 


I\{X, mjj) 


_ \NLOJET+ + 

VNLO[X> m jj 1 _ 

0LO+shower(X, TTljj ) PythiaS ' 


The A-factors decrease with y and thus modify the shape 
of the angular distributions; the impact ranges from a 
few percent at low m 3] to approximately 15% for the 
highest Trijj region. Additional processes accounting for 
electroweak (EW) effects not included in Pythia8 (vir¬ 
tual weak boson exchange and Sudakov-type logarithms) 
are included as EW corrections [35]. The effect is most 
pronounced at high rrijj and low \, and the correction 
factors range from unity at low rrijj to 0.98 - 1.12 in the 
highest rrijj region. The EW corrections and the NLO 
A'-factors are applied as a function of \ an d rrijj to the 
Pythia8 prediction. 

Figure [T| shows the distributions of the data as a func¬ 
tion of x■ The distribution in each rrijj region is nor¬ 
malized to unity, as the sensitivity to new phenomena is 
due to the angular distribution rather than normaliza¬ 
tion. The predicted SM distributions are also shown in 
Fig. a and describe the data well. The EW corrections 
substantially improve the agreement of the SM predic¬ 
tion with data at high rrijj , as can be appreciated from 
the comparison of the predictions with and without these 
corrections shown in Fig. [lj 
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FIG. 1. Distributions of the dijet angular variable x f° r 
several regions of dijet invariant mass rrijj. Shown are the 
data (dots), SM predictions (solid lines) with uncertainties, 
and predictions for a benchmark theory of new contact in¬ 
teractions with a left-chiral color-singlet coupling with con¬ 
structive (dotted, 77 LL = —1 and A = 12 TeV) or destructive 
(dashed, 7 ;ll = +1 and A = 8 TeV) interference with the SM 
processes. Also shown is the impact of the electroweak (EW) 
corrections described in the text. The error bars on data rep¬ 
resent the experimental and statistical uncertainties added 
in quadrature, with a tick representing experimental uncer¬ 
tainties only. The theoretical uncertainties are displayed as a 
shaded band around the SM prediction. 


Models of quark compositeness are probed by searching 
for evidence of new interactions between quarks at a large 
characteristic energy scale, A. At energies below this 
scale, the details of the new interaction and potential 
mediating particles can be integrated out to form a four- 
fermion contact interaction model BB described by an 
effective field theory: 


L qq = [?7ll(9l7/x'2 , l)(<?l7 m <?l) (1) 

+ f?RR(9R7/x9R)(gR7 M 9R) 

+ 27R L (gR7M9R)(9L7 Ai 9L)], 

where the quark fields have L and R chiral projections 
and the coefficients r/ ll, ij rr, and ?7rl turn on and off 
various interactions. 

In this Letter, a contact interaction (Cl) model with 
a left-chiral color-singlet coupling (r] ll = ±1) is used as 
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a benchmark model, as many other models of new phe¬ 
nomena have similar predictions for the dijet scattering 
angle \ at large rrijj. Interference of the signal model 
with the SM process qq —> qq is also included. 

Event samples were simulated with both QCD and con¬ 
tact interactions, taking interference into account and us¬ 
ing the same event generator, underlying-event tune and 
PDF as for the SM simulations. Events were generated 
for both constructive and destructive interference with 
A = 7 TeV and A = 10 TeV. The A = 7 TeV sample 
is then used for extrapolation to other values of A, us¬ 
ing the fact that the interference term is proportional 
to 1/A 2 and the pure Cl cross section is proportional to 
1/A 4 . This procedure is validated with the A = 10 TeV 
sample. As with the QCD prediction, a AT-factor cor¬ 
rection is computed to correct the Pythia LO+shower 
prediction to an NLO calculation. Calculations at NLO 
are provided by Cl JET [35] vl.O. 

Uncertainties in the SM and signal predictions include 
theoretical uncertainties and experimental uncertainties 
on the measured angular distributions. Theoretical un¬ 
certainties in the SM and signal predictions are due to the 
choice of PDF, renormalization and factorization scales, 
choice of event generator as well as statistical uncertain¬ 
ties due to limited simulation sample sizes. The impact 
of the uncertainty in the PDF is estimated using Nlo- 
JET++ with three different PDF: CT10, MSTW2008 [36] 
and NNPDF23 37] ■ These uncertainties are negligible 
(< 1%), as the choice of PDF largely impacts the total 
cross section rather than the angular distributions. The 
uncertainty due to the choice of renormalization and fac¬ 
torization scales were estimated by varying those inde¬ 
pendently up and down by a factor two in Nlojet-|— h 
The resulting uncertainty varies with rrijj and x> rising to 
4% at the smallest x values at high rrijj . The uncertainty 
due to the choice of generator is estimated by comparing 
the predictions from the NLO generator Powheg |33j 
vl.O with those of Pythia8 with AT-factors applied. The 
largest uncertainty due to choice of generator is at the 
lowest rrijj values, where it approaches 20%, while for the 
highest rrijj values and smallest x, it ranges from 10% to 
14%. The uncertainty due to the choice of the shower¬ 
ing model is estimated through comparison of Powheg 
samples showered and hadronized with Pythia8 v8.175 
to Herwig [32] v6.520.2 samples using Jimmy [JO] @T| 
v4.31. The largest value of this uncertainty is less than 
1% at the highest rrijj values and smallest values of X- 
Finally the statistical uncertainties on the A'-factors due 
to limited simulation sample size are small and set to 1%. 

The experimental uncertainty is dominated by the in¬ 
dependence of the jet energy scale calibration. This un¬ 
certainty varies from approximately 15% at small values 
of x for the highest rrijj values, to a few percent at lower 
rrijj values and higher \ values. The uncertainty in the 
beam energy is found to introduce a negligible contribu¬ 
tion. The total uncertainty at the lowest x, highest rrijj 
amounts to 20%, decreasing to a few percent at high \- 
The total theoretical and experimental uncertainties are 


shown in Fig. [T] 

The p-value for the SM hypothesis is (0.25) 0.30 for 
the (second) highest rrijj bin. In the absence of signifi¬ 
cant deviations from the SM prediction, upper bounds on 
Cl contributions are calculated using a one-sided profile 
likelihood ratio and the CLg technique @2] @2], evaluated 
using the asymptotic approximation |44j on events with 
rrijj > 3.2 TeV; the validity of asymptotic approxima¬ 
tion was confirmed using toy simulations. These bounds 
exclude a compositeness scale below 8.1 TeV in a de¬ 
structive interference scenario and below 12.0 TeV in a 
constructive interference scenario. The median expected 
limits are 8.9 (14.1) TeV for the destructive (construc¬ 
tive) interference scenario. 

In summary, dijet angular distributions have been 
measured by the ATLAS experiment in 17.3 fb _1 of 
8 TeV pp collisions at the LHC. Over a wide angular 
range and dijet invariant mass spectrum, the data are 
well described by QCD predictions at NLO. A model of 
quark compositeness is used as a benchmark for theories 
of new phenomena that include new forces and mediating 
particles; such theories predict deviations at small values 
of x- A compositeness scale below 8.1 (12.0) TeV in 
a destructive (constructive) interference scenario is ex¬ 
cluded at 95% confidence level, similar to results from 
the CMS Collaboration m and representing a signifi¬ 
cant enhancement in sensitivity relative to the previous 
limit (at 7.6 TeV for destructive interference) from the 
ATLAS Collaboration |TT] . 
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